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Abstract 

The critical current density of a superconducting film can be easily determined by an inductive 
and contactless method. Driving a sinusoidal current in a single coil placed in front of a super- 
conducting sample, a non zero third harmonic voltage V3 is induced in it when the sample goes 
beyond the Bean critical state. The onset of V3 marks the value of current beyond which the sample 
response to the magnetic induction is no more linear. To take into account, in a realistic way, the 
magnetic coupling between the film and the coil we have developed a discrete model of the inducing 
and induced currents distribution. In the framework of this model the magnetic field profile on the 
sample surface and the coefficient linking the current flowing in the coil and the critical current 
density Jq of superconducting thin films is evaluated. The numerical results are checked measuring 
Jc of several thin films of Y Ba2Cuz0 7 _ s of known superconducting properties, used as a control 
material. 

PACS numbers: 74.25.-q,74.25.Fy,74.25.sv,74.25.Nf 
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I. INTRODUCTION 



The use of superconductors in engineering applications like high field superconducting 
magnets is of fundamental importance. For this reason the critical current density is a rele- 
vant figure of merit, that needs to be measured in samples of different nature and shapes like 
sintered pellets, tapes or thin films. In particular, the development of flat cables for super- 
conducting magnets requires a flexible and reliable measurement of this transport property 
without damaging the sample. 

Years ago Claassen et al.— proposed a non destructive and contactless method to measure 
Jc in superconducting thin films. This method was theoretically justified^, subsequently 
extended for the analysis of bulk superconductors^, and then used for the determination 
of the current-voltage characteristics of superconducting thin films^. The basic principle is 
the following: a small coil, placed in close contact with the sample under test, is excited 
injecting a sinusoidal current Iocosut, and an ac magnetic field is generated. The linear 
and nonlinear response of the superconducting film is then collected by the same coil. The 
critical current density Jc is determined by looking at the behaviour of the induced third 
harmonic signal VscosiSuit) as a function of the inducing current. The link between Jc and 
the current circulating in the probe coil is found^ by assuming two important conditions: i) 
the sample surface is flat and extended to infinity in two dimensions; it) the total induction 
field is zero on the back side of the sample. Of course this last condition follows from the 
first one when the magnetic penetration depth A < d, where d is the sample thickness. 
Due to these assumptions, the field profile must have as boundary condition that only 
the parallel component is nonzero at the sample surface. By using a "mirror technique", 
the sample is replaced by an image coil reflected by the sample surface and carrying the 
same current. In this case the total induction field on the sample surface is just twice the 
value that would be generated by the single coil in the absence of the screening currents. 
According to %%) and the Ampere's law, the maximum sheet current K s that can be sustained 
by a superconducting film placed at a distance h from the coil results to be Ksih) = 
j^Bj OTAL (h) = -^2 ■ max{B T {h)} where B r {h) = fi I f(r, h)cos(ut) is the radial component 
of the induction field generated by the coil at distance h in the absence of the sample. 
Applying the Bean's model, Ks = J c d. By equaling the two expressions found for Ks, the 
scaling relation between the critical current density and the inducing coil current follows: 
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J< = k.!f (!) 

where Icq is the inducing critical current (defined, for example, using a threshold criterion) 
and k = 2 * max{f(r, h)} is the coil factor. 

In this paper we review this contactless technique by using a discrete modelisation for 
the current distribution in the coil, in order to calculate in an exact way the induction field 
in all the space as a function of the sample distance to the coil, thickness and magnetic 
penetration depth. We verify how and under which circumstance the assumptions made are 
verified in the framework of the model, and then we calculate the scaling coefficients linking 
the current circulating in the coil to the critical current density of the sample. The paper is 
organized as follows: in section II, after a brief overview on previous studies, we present the 
model and calculate the vector potential in all the space for the sample coil configuration; 
in section III we check the validity of our model comparing our analysis with similar results 
presented in literature; in section IV we present the experimental set-up used to measure 
superconducting samples of known transport properties to validate the analytical results; in 
section V we finally draw some conclusions. 

II. THE DISCRETE MODEL 

The Non Destructive Evaluation (NDE) of the conducting properties of metallic spec- 
imens by using eddy currents methods has been since years an important subject for the 
scientific community for both applicative and fundamental issues. As an example, we can 
cite the pioneering work by C. P. Bear£ that used a simple solenoid to measure the resistivity 
of a sample placed inside. The necessity of studying large area samples raised the problem 
of calculating the spatial profile of the induction field when using a small multiturn coil. 
This problem was solved in an original way by Dodd and Deeda^ with an "hybrid" model: 
they calculated the field profiles for a single turn coil and extended the result to the case of 
multiturn coils by integrating over the coil section. Their work has been the starting point 
for hundreds of experimental and theoretical studies in the field of NDE. 

In the case of superconducting materials, the thin film technology raised the interest in the 
study of both the transport properties^ (the critical current density) and the electrodynamic 
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response (the absolute magnetic penetration depth) 8 i 9 i 10 i n i 12 i 13 i 14 i 15 without damaging the 
sample under test. For these issues, a two coil configuration with both continuous and 
discrete approach has been used. 

For a single coil configuration, a completely discrete approach by Gauzzi et ai^ was used 
to model the coil impedance variation for the measure of the magnetic penetration depth 
of superconducting samples of different nature and shapes. Recently M. W. Coffey^ has 
reviewed the single and multiturn impedance coil problem by using a continuous method of 
kernel functions. For the evaluation of the transport properties of superconducting samples, 
we have shortly presented in the introduction the continuous model approach^. 

In this section, we extend the discrete model by Gauzzi et ci1f& to the calculation of the 
field induction profiles in all the space. We will show how and in which cases this exact 
model is simpler and more precise respect to previous hybrid or continuous models cited 
above. 

We investigate the ac response of a superconducting slab of finite thickness situated at 
— d < z < 0, at a distance h from a pancake coil, as shown in figure [TJ Further, we consider 
the hypothesis that the outer radius of the coil is less than half sample size to neglect edge 
effects. This will be better clarified in the following, when we will evaluate and show the 
magnetic field profiles. 

First, we model the coil as a discrete distribution of concentric turns axially and radially^ 
equispaced. The wire is supposed infinitely thin; this assumption is valid when the wire 
diameter is much smaller than the skin depth at the working frequency. For example, in 
the case of a copper wire and at frequencies below MHz this approximation holds for wire 
diameter up to 0.1 mm. Thus the total current density circulating in the coil can be written 
in cylindrical coordinates in the following way: 

N-l M-1 

Je,coii(r, z) = I^2$(z + h + nAh) ^ 6(r - R - mAR) (2) 

n=0 m=0 

where I is the total current flowing in the coil, S is the delta function, R is the internal 
coil radius and Ah, AR are the spacings between N adjacent sheets and M adjacent turns 
in the same sheet respectively (see fig.l). We assume that Ah and AR both equals the wire 
diameter d w . 

Secondly, we consider the Maxwell equation at a frequency ui /2tt for the vector potential 
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A in the London gauge (V ■ A = 0), explicitly separating the contributions to the current 
density in the coil and in the sample: 



V A^j — fJ>o(juj,coil ~\~ j ui , sample) ■ (3) 

At low frequencies within the limit of validity of Ohm's law, the induced screening current 
in the sample can be expressed in the following way: 



jui, sample i^O' u) A UJ (4) 

where a u is the frequency dependent complex conductivity of the sample. By inserting 
eq. lUinto eq. C3 the solution for the field A u becomes a function of j^^ou and a^, the latter 
being the only free parameter. 

To solve eq. El we note that the current density ] co u in the coil lies in the (x-y) plane 
of the sample. The same must be true for the induced screening currents in the sample, 
as in a mirror image. This consideration, and the assumption of a sample with isotropic 
conductivity a in the plane, allow us to express eq. Elm cylindrical coordinates. In this case, 
only the orthoradial component A# of the vector potential A is non zero: 



V 2 A? = -fj, (ju t coii + iH d (z)u<T u Ao) (5) 

where the origin and direction of the z-axis are chosen as shown in figure [U the u- 
subscripts have been omitted, and H d (z) is the Heaviside function, whose value is unitary 
in the [0,d] interval and zero elsewhere. We have solved the differential equation El in all 
the space with and without the sample by following the Pearl's formalis m 1 ^ . The solution 
for the vector potential without the sample is: 



A#(r,z) = -A*oy J 



|z+/l n |7 



^]f m Ji(r m 7) 



Ji(r7)c?7 



(6) 



In the presence of a superconducting sample we have three different solutions, corre- 
sponding to three different regions: 



z<0 



Aa(r,z) = /i y J 



;i-F( 7j ^d,0) Ji(ry)d7 (7) 



1 + 2 7 2 / 2 + 27/^/1 + / 2 7 2 coth + Z 2 7 2 ) 



< z < d 



(8) 



A^(r, z) = — /x 



2 Jo 



y^exp(-/i w 7) y^r m Ji(r m .m7) 



G( T ,z,d,Z)J 1 (r 7 )d7 (9) 



z>d 



27/ 


7/ sinh 




-y/l + I 2 J 2 COsh (^fV 1 + / 2 7 2 y 




(1 + 2 7 2 / 2 ) sinh + Z 2 7 2 


) + 2 7 Zyi + Z 2 7 2 cosh Uyjl + l 2 ^ 


2 ) 



(10) 



A#(r,z) 



<o 



T r®o 
if] 



[1 



-Z7 



L( 7 ,Z,d) 



2 7 /v/l + ZV 



L(>y,l,d)e-*r] JiM^ll) 



(12) 



;i + 2 7 2 / 2 ) sinh (f^l + / 2 7 2 ) + 2 7 / v /rT7 2 7 2 cosh (fyi + Z 2 T 2 ) 

where 7 is a wave number in the xy-plane, the functions F, G and L represent the sample 
contribution in the considered region, and Jo and J\ are the Bessel functions of order zero 
and one respectively. I is a complex length expressed by the formula: 



/ 



(13) 



Note that it reduces to the normal skin depth for the case of normal metals and to the 
magnetic penetration depth for the case of superconductors. The radial and axial com- 
ponents of the magnetic field with and without the sample have been calculated from the 
vector potential by using the standard recursive relations: 
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(14) 



B z {r iZ ) = ~[rM) (15) 
r or 

The explicit expressions for the components of the induction field in all the space are 
reported in appendix A. 

III. NUMERICAL RESULTS 

Our general model has been checked in the particular case of a single turn, single sheet 
coil by using the calculations of Dodd and Deeds^. In figure El we compare the results of the 
calculations of the radial (B r ) and axial (B z ) components of the induction field in the case 
of no sample for a single turn coil with inner radius r =l mm and wire diameter equal to 50 
fim. The fields have been calculated at a distance z = r from the coil surface (z — 0). The 
two models give the same result. In case of multi-turns coils, the continuous approach^ 
has been widely used. We point out that assuming the simple rectangular coil section of a 
continuous model produces an over-estimation of the effective area covered by the exciting 
coil current in respect to the discrete case, in which the effective surface is the sum of N-turn 
discs of the diameter of the coil wire. One would expect therefore an under-estimation of the 
maximum value of the magnetic field compared to the discrete case; the error decreases as 
the number of turns increases approaching the continuous limit. To show this, we compare 
the continuous^ and the discrete model calculating the radial component of the magnetic 
field in two cases: a) a coil of 5 turns, 5 sheets; b) a coil of 22 turns, 24 sheets. The inner 
radius ro is 1 mm and the sample to coil distance is h=0.1 mm; the wire diameter has been 
assumed to be 50/im. The results of these calculations are reported in figure 03 As expected, 
the maximum B r field calculated with the continuous model is under-estimated respect to 
the discrete one of 8 and 2% for the case a) and b) respectively. By using the methods 
proposed in ref;^, this would imply an under-estimation of the critical current density of 
the same amount. 

Once the discrete model has been successfully tested, we have used it to calculate the 
coefficient k(h) in eq. [TJ Firstly, we have verified the condition to apply the Bean model 
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properly. To do so, we have calculated the total induction field on the sample surface by 
using eq. |H1 and HU The parameters used in this calculation are: the sample thickness 
d=700 nm, its magnetic penetration depth A = 250 nm for a standard high quality YBCO 
sampl e 1 '^ 1 '? at 77 K. For the probe coil we have considered a 22 turns, 24 sheets coil with the 
inner radius r = 1 mm and a wire diameter equal to 50 /im. The results of this calculation 
are reported in figure the normal component B 2 is negligible in respect to the radial 
component B r . Moreover, it is worthwhile to emphasize that the magnitude of B r is about 
two times larger than in the case without sample shown in figure 02 in the same point of 
the space, as we can expect within the framework of a mirror image coil solution analysis. 
For the Bean model be valid, the total induction field must be zero on the opposite side 
of the film. This is true when the effective screening length \ SC reenin g is equal to or less 
than the film thickness d. From eq. 02 in the case of a superconducting sample, it results 
^screening = A • coth(d/ A) 16 i 2Q . When X S creening > d the error in the estimation of Jc may 
be very large: to estimate it roughly, we can calculate the relative difference in the Jc 
values when the film thickness equals d and the effective screening length X S creening, i-e. 
|1 — Jc (^screening) / Jc(d)\- F° r example, for a film with <i=200nm and A=250nm, the error 
can rise up to 46%. 

To calculate the coefficient k we take the maximum of the coil function f(r) plotted in 
figure 03 This represents the spatial region where the field reaches its maximum value. In the 
case of the copper coil reported above, our discrete model gives the value of k = 2 * 10 5 m _1 
at a sample-coil distance h=0.1 mm. Our calculation allows to estimate k(h) with an error 
of 2% for a minimum 10% error in the determination of h. k(h) decreases as h increases 
as it is shown in figure El because of the presence of terms oc e~ hl in the field expression. 
In the same graph, the experimentally determinated values of k at various distance are 
also reported. The experimental values are in good agreement within the error bars with 
the behavior predicted by the discrete model. Details of the measurement will be given in 
the next section. The exact knowledge of the relation linking k to h makes this technique 
extremely flexible. For example, in the case of superconducting tapes or large area films, 
the probing coil must be moved all along the sample. Thus, to correctly estimate Jc, it is 
necessary to know exactly the magnetic coupling versus the sample to coil distance h. 
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IV. EXPERIMENTAL 



To validate our discrete model we have inductively measured the critical current density 
of a superconducting film of well known transport properties. For the measure of the third 
harmonic signal V3 we used a Lock-in amplifier (Signal Recovery 7265). The amplitude 
of the inducing sinusoidal current (the Lock-in reference signal) circulating in the coil is 
measured by using a high precision resistor in series with the coil. A commercial highly 
linear amplifier increases the inducing sinusoidal signal. The output signal is filtered by 
using a home made, low noise notch filter to reduce the first harmonic component noise and 
to avoid any saturation of the Lock-in. The temperature is monitored by a temperature 
controller (Neocera2021) connected to a Cernox resistor thermometer. The characteristics 
of the coil used in our experiment have been reported in the previous paragraph. 

The samples tested are 700 nm thick YBCO thin films produced by Theva GmbH with 
dimensions 10 x 10 mm 2 . The critical temperature and the critical current density at liquid 
nitrogen temperature given by the producer are 87.9 K and 2.57 MA/cm 2 (expressed in 
root mean square) respectively. The value of the critical current Icq was experimentally 
determined looking at the onset of non linearity by the extrapolation of the linear portion in 
the V3 versus I plot (see for example in figure M and ref.— ). The value found for the critical 
current density of the sample of fig. H3 is Jc = (2.6 ± 0.3)MA/cm 2 at 77 K, in very good 
agreement with the data reported by the producer. For all the other samples tested, Jc has 
been found statistically well within the error bars of the measurement presented above. The 
error on Icq is estimated as the difference between the first point that is no more aligned 
on the linear portion V3 versus I and the extrapolated Icq value (see the arrows on figure 
EJ). This error is larger at low temperature because of the increase of the critical current 
density: at higher Jc, the noise level introduced by the amplifier is enhanced thus producing 
a rounding in the lossy state transition. 

We have repeated the same experiment at different sample-coil spacings h, by using a 
micrometric screw with 10 fim resolution placed under the sample holder, to calculate the 
experimental scaling factor k vs h and comparing it to the expected theoretical behaviour. 
We note some important points: i) the experimental values have been obtained by consider- 
ing the critical current density equal to the value given by the producer; it) to minimize the 
error in the determination of the sample-coil distance, we have set the zero of the micromet- 
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ric screw when the sample is in contact with the coil (h = in) the error bars for k(h) 

have been obtained by taking into account all the possible error sources (d and Icq) whereas 
for h we have taken the minimum resolution of the micrometric screw. As one can see in 
figure El there is a very good agreement between the k(h) values predicted by the discrete 
model and the experimental ones. 

As a final test, we have plotted Jc versus T close to T c in figure The behaviour of the 
critical current density is a linear function of the temperature. This is usually explained as 
due to the presence of Josephson coupling between grain o 21 ! 22 . A more detailed analysis of 
the Jc(T) behaviour for T > Tc/2 goes beyond the scope of this work and will be presented 
elsewhere. 



V. CONCLUSIONS 



In this paper we have revisited in detail an inductive and contactless technique used for 
the measurement of the critical current density in superconducting films. The inductive 
geometric coupling factor has been calculated exactly by developing a precise and realistic 
model describing the current circulating in a real coil, and its image in a superconducting film 
placed in close contact to it. In particular, the scaling law of the geometric factor versus 
the sample-coil distance is found. This is a fundamental issue to render this technique 
flexible, in order to make possible a local non destructive evaluation of the critical current 
density in large areas films or superconducting tapes just moving a single probe coil at 
a known distance from the testing sample. This discrete analysis has been successfully 
verified by a comparison with similar continuous models presented in literature. The model 
presented here is more realistic, relatively simple and can describe the experimental data 
with a higher degree of precision in particular in the case of coils with a low number of 
turns contrarily to what happens using a continuous approach. The model has been verified 
experimentally by measuring the critical current density on YBCO films having well known 
transport properties. 
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APPENDIX A: THE INDUCTION FIELD EXPRESSIONS 

By using expressions ( l6llT2| , ( fTIlfT5|) and the recursive relation for the derivative of the 
Bessel functions we have calculated the radial and the normal induction field in all the space 
with and without the sample. The expressions are the following: 

without the sample: 



B r (r,z) = /x y J 



^ e -|^l7^ rm J l(rm7 ) 



Ji (7-7)7^7 



(Al) 



B z (r, z) = nJ± 



^ e -k+M7^ rmJi(rm7) 



One can define a geometric coil function^ f(r) for the radial component: 



(A2) 



m 

with the sample: 

z<0 



Ji ^7)7^7. 



(A3) 



B r (r, z) 



t r°° 



^ e H^l7^ rm j l(rm7) 



(l-F( T ,d,0) Jiinhd-y (A4) 



B z (r,z) = /i : 



(1-F( 7 ,d,0) J (r 7 ) 7 d7 (A5) 
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< z < d 



B r (r, z) = ^ 



n m 



<9G( 7 , Z, d, I) 

dz 



Ji(n)d-f (A6) 



BJr, z) 



n m 



G(-f,z,d,l)J (n)-fd-f (A7) 



dG(>y,z,d,l) 27 ^^T+ZVcosh (^^MV) - (1 + /V)sinh (^v 7 ^^ 



dz 
z>d 



1 + 2 7 2 / 2 ) sinh (fv/l + ZV) + 2 7 /v/l + / 2 7 2 cosh (f^l + l 2 r 



B r (r,z) = /i y J 



n m 



e 27 J^r^jdj 



(A9) 



£ 2 (r, z) = -/^oy y 



n m 



[1 



-«7 



L(7,/,rf)e- d7 ] J (r 7 ) 7 #10) 



where F( 7 , z,d,l),G(j, z,d,l) and L(j,d,l) have been defined in eq. [HI eq. [THland eg. IT2l 
respectively. 

Finally, it is possible to evaluate exactly the maximum electric field induced on the surface 
of the sample (z=0). By using the third Maxwell equation and eq. El we have: 



Etf(r, z) = max 



9A$ 
' dt 



fioloujmax {f(r)} 
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FIG. 1: Section of the sample-coil configuration. The coil is modelled as a discrete distribution of 
concentric turns axially and radially equispaced. h is the sample-coil distance and d is the sample 
thickness. Ah and AR are the spacings between N adjacent sheets and M adjacent turns in the 
same sheet respectively. 
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FIG. 2: (Color Online) a) radial component B r (r,z = ro) of the induction field for a one turn one 
sheet coil calculated by using the Dodd and Deeds's model (continuous line) and our discrete model 
(open circles) respectively as a function of r/ro; b) the same as in a) for the normal component 
B z (r, z = r ). 
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FIG. 3: Radial component of the magnetic field for the two coils of cases a) and b). The open 
circles and the lines represent the discrete and the continuous model respectively. 
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FIG. 4: Radial component of the total induction field (open squares) for a multi-turn coil. In the 
inset the result for the normal component (open circles) is shown. Both components have been 
calculated at the sample surface as a function of t/tq. 
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FIG. 5: (Color Online) Determination of the scaling factor k versus the sample to coil distance 
h. The continuous and the dashed lines represent the behaviour expected for the discrete and the 
continuous model respectively. The open circles represent the fc-values experimentally evaluated. 
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FIG. 6: Determination of the critical current Icq- The arrow indicates the last point below which 
the third harmonic voltage is no more linear as a function of I. 
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FIG. 7: (Color Online) Critical current density Jc* versus temperature for a YBCO thin film. 
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